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We study the possible charmonium-like states with J PC = __ ,0 _+ using the tetraquark inter- 
polating currents with the QCD sum rules approach. The extracted masses are around 4.5 GeV 
for the __ charmonium-like state and 4.6 GeV for the 0~ + charmonium-like state while their 
bottomonium-like analogues lie around 10.6 GeV. We also discuss the possible decay, production, 
and experimental search of the 0~~ charmonium-like state. 

7— I ' 

PACS numbers: 12.38.Lg, 11.40.-q, 12.39.Mk 
i Keywords: Charmonium-like states, QCD sum rule 

I—* 

1 ; I. INTRODUCTION 

'"^ ' Since the Belle Collaboration observed the narrow state X(3872) on the threshold of D°D° in the B+ -)■ 
II ■ K + X(3872) — > K + J/ipTT + Tr~ channel in 2003 [if, many new charmonium or charmonium-like states have been 
Oni observed such as F(3940), F(4260), £(3930), X(3940), y(4325), y (4360V y(4660), £+(4430), Z+(4050), Z+(4250) and 
y(4140) For experimental reviews, one can consult Refs. [12T - fl5l ]. 

The discovery of these new states have enriched the charmonium spectroscopy greatly. It is very difficult to 
accommodate all these states in the conventional quark model. In order to study their underlying structure, many 
interpretations were proposed such as the hybrid mesons, the molecular or tetraquark states, baryonium states and 
so on. For example, X(3872) was speculated to be a hybrid charmonium state in Ref. [l6|, a DD* molecular state 
in Ref.[l7| and a cqcq tetraquark state in Ref. [l8|. Z + (4430) was assigned as a D*Di molecular state in Ref. [l9| 
and a cscs tetraquark state in Ref. [20| . V(4260) was proposed as a hybrid charmonium state in Ref. [2l[ and a 
. cqcq tetraquark state in Ref. [22[- F(4140) was proposed as a D^D* molecular state in Ref. [23|. However, one 
should be very cautious that the conventional charmonium spectrum may be distorted if one considers either the 
coupled-channel effect or the screened linear confinement force [24| . More charmonium states can be accommodated 
■ below 5 GeV within this picture [24| . 

The possible ccqq states with various quantum numbers including J PC = were investigated extensively in Ref. 
[25| in the early 1980s. The author discussed the spectroscopy, decay and production of the ccqq systems with the 
angular momentum L > 1 by taking account of the color magnetic and electric forces. In the conventional quark 
model, states with J PC = arc exotic states. They cannot be composed of a pair of quark and antiquark. In Ref. 
[26( | . we noticed that the light tetraquark currents with J PC = do not support a low-lying resonant signal below 
2 GeV. Since increasing the quark mass reduces the kinetic energy and thus may help stabilize the system, we will 
study whether the ccqq states with J PC = exist in this paper. 

The paper is organized as follows. We construct the tetraquark currents with J PC = 0~~,0 ^ using the diquark 
and antidiquark fields in Sec. II and derive the spectral densities in Sec. Ill which are collected in the Appendix. 
We perform the numerical analysis and extract the masses in Sec. IV and discuss the possible decay, production, and 
experimental search of the charmonium-like states in the last section. 

II. TETRAQUARK INTERPOLATING CURRENTS 

We have constructed the light tetraquark interpolating currents with J PC = using the diquark- antidiquark 
fields in the previous work (2(|. In this work we follow the same steps. We first construct ten color singlet pseudoscalar 
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operators considering both the Lorentz and color structures: 56, P§, Vq, Aq, Tq, S 3 , P3, V3, A 3 , T 3 . The subscripts 6 and 
3 indicate that the color structures of the tetraquark are 6 (£> 6 and 3 Cg> 3 respectively. Details can be found in Ref. (2(| . 
With the charge-conjugation transformation we get 

C56CT 1 = V 6 , CA 6 C _1 = P 6 , = P 3 , C5 3 Cr 1 = V 3 , CTeCT 1 = T 6 , Cr 3 C _1 = T 3 . (1) 

Using the above charge-conjugation relations, we can obtain the currents with definite C parity: 

1. For the quantum number J PC = : 

T)i = 5 6 - V 6 = q^Cctiqa^Ccl + qblhCcl) - q 7 a ' C-f 5 c b (q a Ccl + q b Ccl) , 

n 2 = A 6 -P 6 = qaC-f^CbiqaY^Ccl + q b ^ 5 Cc^) - qlC-i^c b {q a YCc{ + q b ^Cc T a ) , 

m = A 3 -P 3 = q^C^c b (q a ^ l5 Ccl - q b ^Cc T a ) - q T a C llilb C b {q a ^Ccl - %^Cc T a ) , (2) 

V4 = S3 - V 3 = qlCc^qa^Ccl - qhl^Ccl) - q\ ' C"f 5 c b (q a Ccl - (ftCcf ) . 

2. For the quantum number J PC = h 

V5 = 5 6 + V e = qaCc b (q a 'y 5 Cc b r + qblsCc 7 ) + q 7 C^ 5 c b (q a Ccl + q b Cc 7 ) , 

Ve = T 3 = qaCa^Chiqaa^-fsCcl - qhcr^^Cc 7 ) , 

T) 7 = A 6 + P 6 = qlC lllCb {q a ^Ccl + q^l^Ccl) + qlC llil5 c b {q al »Ccl + q^Cf?) , 

V8 = A 3 +P 3 = qJC^c b (q a n 5 Ccl - q bl ^Cc T a ) + q T a C W C b {q a ^Ccl - q b ^Cc T a ) , (3) 

% = S 3 + V 3 = qaCc b (q a -f 5 Ccl - q^Ccl) + q J a 'Cj 5 c b (q a Ccl - q b Cc 7 ) , 

VlO = T 6 = qlCo^CbiqaO^^Ccl + q b o^ l5 Cc T a ) . 

It is understood that Eqs. ©-© should contain (ucuc + dcdc) in order to have definite isospin and G-parity. Because 
of the SU(2) flavor symmetry, we do not differentiate the up and down quarks in our analysis and denote them by q. 



III. QCD SUM RULE 

QCD sum rule is a powerful approach to study the hadron properties in the past several decades [271429I ] . We consider 
the two-point correlation function: 

n^ 2 )^ [ d A xe iqx (0\Tri(x)rf(0)\0), (4) 



where 77 is an interpolating current. At the hadron level, the correlation function Tl(q 2 ) is expressed via the dispersion 
relation: 

/■oo n phen( \ 
W hen [p 2 )= / P ISj 

Jo s-p z - is 

The spectral function reads: 

pPhen{s) = ^ (5(s _ TO 2 )(()|??|7l)(?l|r;+|0) 

n 

= fxS(s — m 2 x ) + continuum, (6) 

where the usual pole plus continuum parametrization of the hadronic spectral density is adopted, fx is the overlapping 
parameter of the current to the pseudoscalar state X: (0\n\X) = fx- 

At the quark- gluon level, the spectral density can be evaluated up to dimension eight with the same method in 
Refs.[l8|, Ho, l30l 1 3 1| - Omitting the light quark mass, we use the coordinate-space expression for the light quark 
propagator and momentum-space expression for the charm quark propagator: 

AXab ■ \n 1 sab rah 2 

^ = ^ + 32^f "-^K* + **»») - -(qq) + —(q 9s O- ■ Gq), 
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where x = lfl x^ p = (qg s a ■ Gq) = (g s qa^G^q), (g 2 GG) = ( 3s 2 G^G^), {glfGGG) = (g 3 J abc G a 7S GlG^), a 

and b are the color indices. The momentum-space propagator with three soft gluon lines can be found in Ref. [28j . 
For the light quark propagator, we use the D dimension coordinate-space expression. The dimensional regularization 
is used throughout our calculation. The Tl(q 2 ) in the operator product expansion(OPE) side can be written as: 



n 



OPE 



(P 2 ) 



X>PE 
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Performing Borcl transformation for the correlation function, we arrive at: 

/| e -™*/Mi = H d se- s ^p OPE (s). 

where sq is the threshold parameter. The mass Mx reads: 



(8) 



(9) 



m x = 



■so i - 
4m? dse 



s/Ml p OPE^ 



(10) 



For all the tetraquark currents in Eqs. (|2|) and ([3]), we collect the p (s) in the Appendix. Both the quark 
condensate (qq) and quark gluon mixed condensate {qg s <? • Gq) vanish in the chiral limit m q = 0. One may wonder 
whether the quark condensates proportional to the charm quark mass: m c (qq) and m r (qq ^a ■ Gq) exist. They are 
usually very important corrections in the scalar, vector and axial- vector channels [l8l. l3l[ |3~I1. |35|. However, these 
terms also vanish in the pseudoscalar channel. The diagrams in FigJT] vanish due to the special Lorenz structures of 
the currents. 





FIG. 1: Feynman diagrams for m c {qq) and m c {qg a a ■ Gq). In each diagram, the two upper lines represent the light quark 
propagators and the two lower lines represent the charm quark propagators. 
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IV. NUMERICAL ANALYSIS 



We use the following values of the quark masses and various condensates in the QCD sum rule analysis [14112 
The charm and bottom quark masses are the running masses in the MS scheme: 

m c (m c ) = (1.23 ±0.09) GcV, 
m b (m b ) = (4.20 ± 0.07) GcV , 
(qq) = -(0.23±0.03) 3 GcV 3 , 
(qg s <r ■ Gq) = -M^qq) , 
M 2 = (0.8 ± 0.2) GeV 2 , 
(ss)/(qq) = 0.8 ±0.2, 
(g 2 GG) = 0.88 GeV 4 , 
(glfGGG) = 0.045 GeV 6 . 



(11) 
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FIG. 2: The variation of Mx with Mg (left) and so (right) for the current rj2. 



The Borel mass Mb and the threshold value so are two pivotal parameters in the numerical analysis. The working 
region of the Borcl parameter is determined by the convergence of the operator product expansion and the pole 
contribution. The requirement of the convergence of OPE leads to the lower bound Mf nin of the Borel parameter. 
The upper bound M^ ax of the Borel parameter results from the requirement of the pole contribution. 

The four quark condensate (qq) 2 is negative and the most important condensate correction numerically. Its absolute 
value is much bigger than the gluon condensate (g 2 G 2 ) and the dimension 8 condensate (qga ■ Gq)(qq). Technically, 
we require that the four quark condensate be less than one third of the perturbative term to ensure the convergence 
of OPE, which results at the lower bound of the Borel working window, M? nin ~ 2.4 GeV. 

The pole contribution (PC) is defined as 



(12) 



which depends on both the Borel mass Mb and the threshold value Sq- sq is chosen around the region where the 
variation of mx with M B is the minimum. For example, we choose sq ~ 25GeV 2 from the variation of the mass with 
so as shown in Figs. OandO Requiring that the PC is larger than 40%, one gets the upper bound M^ ax of the Borcl 
parameter M B . We list the working range of the Borel parameter for the ten tetraquark currents 771 ~ 7710 in table HI 
The masses are extracted using these threshold values and M B = 3.5GeV 2 , which are also listed in Table HI Only the 
errors which arise from the uncertainty of the threshold values and variation of the Borcl parameter are taken into 
account. Other possible error sources include the truncation of the OPE series and the uncertainty of the condensate 
values etc. The last column is the pole contribution. 
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Currents 
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4.67 ±0.10 


56.8 



TABLE I: The threshold value, Borel window, mass and pole contribution for rji — 7710. The mass and pole contribution are 



calculated at M f 



3.5GeV" 



From the variation of the mass with sq, there is a plateau in the region of sq = 9— 13GeV 2 for r]2, 774, 776, 779, Vw as 
can be seen in Fig. [2] This plateau looks like a resonance signal. However, it's just an unphysical artifact because 
both the numerator and denominator in Eq. (|10p are negative within this region. The variation of Mx with the Borcl 
parameter is weak. 

For 771, 773, 775, 777, ry 8 , the extracted mass Mx grows monotonically with the threshold value. Also the variation of 
Mx with the Borel parameter is significant. So we do not present the numerical values in Tabic Q] for these currents. 
These currents may couple to the ,0 states very weakly. The continuum contribution may be quite large. 
These two factors may lead to the above unstable mass sum rules. 
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FIG. 3: The variation of Mx with M% (left) and so (right) for the current ?76. 

Replacing m c with mi, and repeating the same analysis procedures done above, wc collect the results of the 
bottomonium-likc systems in Table ITT1 
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TABLE II: The threshold value, Borel window, mass and pole contribution for r/u — 77106- The mass and pole contribution are 
calculated at Mg = 9.0GeV 2 . 
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V. SUMMARY 

We have constructed the charmonium-like tetraquark interpolating currents with J = ,0 h using the 
diquark-antidiquark fields. Then we calculated the correlation functions and the spectral densities of these currents. 
Both the dimension 3 quark condensate and dimension 5 quark gluon mixed condensate vanish if we take m u ( j = 0. 
The special Lorenz structures of the currents prohibit their appearance in the OPE. The four quark condensate (qq) 2 
becomes the most important power correction numerically. It is much bigger than the gluon condensates (g 2 GG), 
{g 3 fGGG) and the dimension 8 condensate {qga ■ Gq){qq). 

In the working region of the Borel parameter, the variation of the extracted mass with sq and Mb is stable for the 
currents 772, ??4, r]6i V9i Vw- F° r the charmonium-like states, its mass is around 4.5 GeV. For the ^ charmonium- 
like state, the mass is around 4.6 GeV. For the and h bottomonium-like states, their masses are around 10.6 
GeV. It's interesting to note that the extracted mass value ~ 4.5 GeV of the charmonium-like state is quite close 
to the mass value 4.1 — 4.4 GeV in Ref. [25| . 

The possible decay modes of the charmonium-like state are straightforward after making Fierz transforma- 
tion to the diquark type interpolating currents in Eq. O They can be expressed in terms of the linear combi- 
nation of the meson-meson type of operators such as: (q a l^q a )(cbY l l5Cb), (?a7/i75<7a)(c&7' 1 c&), (q a l^c a )(c b ^j 5 q b ) + 
(q~aJn75C a )(cb^qb), (q a c a )(cbl5qb) - (q a j5C a )(c b q b ) . 

There are two types of the 0~~ charmonium-like state with different isospin and G-parity: I G = 1 + and I G = 0~. 
Considering the conservation of the isospin, G-parity and C parity, we collect the S-wave and P-wave decay modes of 
the possible charmonium-like state in Table Hill 

Clearly the S-wave decay modes are dominant. The S-wave decay products always contain a P-wave and S-wave 
meson pair. Such a decay pattern is also characteristic of the hybrid meson. Although the exotic state can not 
be composed of a cc pair, it can be a cGc hybrid state. In fact, the charmonium-like tetraquark operator and 
the cGc hybrid operator probably couple to the same physical state. 

This interesting CP~ state may be searched for experimentally at facilities such as Super-B factories, PANDE, 
LHC and RHIC in the future. Especially at RHIC, plenty of charm, anti-charm and light quarks are produced 
simultaneously which may be helpful to the formation of the charmonium-like state. 





S-wave 


P-wave 


0- 


D*(2007) u L>i(2420) u + cc, 
£>2(2400)°D°(i865) + cc, 
w(782)xci(lP), Wi(1285) 


JJ U (1865)L>*(2007) U + cc, 
L»*(2007)°D*(2007) , 
J/^r,, J/H, ^S) V , Jfc(15)w, 
i lc {2S)u, h c (lP)a, fcc(lP)/ (980) 


1+ 


£>*(2007) u L>i(2420) u + cc, 
D* (2400)°£>°(1865) + cc, 
p(770) X ci(lP), J/M(1260) 


JJ U (1865)L>*(2007) U + cc, 
L>*(2007)°D*(2007)°, 
Vc(lS)p, rjc(2S) P , /i c (lP)a (980) 

J/l/'TT, J/^7Tl(1400), lfj(2S)TT 



TABLE III: The possible decay modes of the charmonium-like state. 
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Appendix A: THE SPECTRAL DENSITIES 



In this appendix we show the spectral densities of the tetraquark interpolating currents defined in Eqs. (0)-([3]). The 
same subscripts are used to denote the results for the currents T)i,r)2,fl3) r }'Li r l5i r l6iV7)' r l8)' r l9i VW- 

p OPE = pP ert^ + p<9 - 9 ) (s) + p <G 2 )( s) + p mix + p {w) 2 ^ ( A1 ) 

For the expressions below, the integration limits are: 



1 + s/i - 4ml /s _ 1 - V 1 - 4m c/ s 



UbUbJu ^ ' llbl/lb ^ 

R - 1 _ R QW c 

Prnax 1 ^7 Pmin 9 ■ 

as — mf. 

1. The spectral densities of the interpolating currents with the quantum numbers J PC = 0~ 
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For 771: 



PT rt (s) = ["'-da /' Pra "^ (1 _ a _ j9)2[(a + j9)m 2_ 3a ^ ][(a + ^ m 2_ a ^ ] 3 ; 



2 7 tt 6 J amte a 3 i^ min /3 3 
= 0, 



(GG), s (ff?GG) [ ama * da r™* , nr (l - a-/3) 2 m r „, „, 2 „ „ 



^{s) = -^fVT^Ts, (A2) 

jjmte^^x = m 2 c (qg s a-Gq)(qq) f 1 da 2m 2 . | - -^^-Sg^ ^ 

n ( GGG) ., 2 , {glfGGG) f 1 [P*™ 

Mgg(l - a - P )M%) - 2 In 2 - lE [12(a/gM| )2 + 6a/W| (q + ^2 + (q + ^ 
ap 

, 33(a/3Af 2 ) 2 + 12a/?M 2 (a + p > 2 + (a + 0fmj + (1 - a - /?) 2 mg + (a + ^ 



ap 1 



2a 3 



[3a/3A^ + A/^(a + /3)m 2 ]}e <^ M S . 



For 772: 



2 5 7r 6 7 Q _ a 3 i^ m . n /? 3 
p^>( S ) = 0, 

4 >M = 4^/ 7^ *P£ 3^[2(a + «m»-3a^* 



pT x {s) = 0, 



TjTmxiqq) , n T 2 , _ m 2 c {qg s (j ■ Gq){qq) f 1 da 4m 2 - a(1 _1) M 2 
n 2 - 7-9 / T~l7~n72 °/ e * 



(A3) 



6tt 2 y a v aM 2 



Mgg(l -»-j)M|) -2h2-7g f| „ + 6q(3m|(q + 2 + 

ap 

33(a/3M 2 ) 2 + 12a/3M 2 (a + /3)mg + ( Q + /j) 2 m^ (1 - a - /?) 2 mg ro fl)lf2 , fl . 2 , 

+ o 1 3 [2apM B + (a + /3)m c 

H — n/ - R\ 2 _(a+£)m| 

^3^[3a/3Ml + Ml(a + ^K]}e "^F . 
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For 773: 



pT r \s) 



pf«\s) = 0, 

jgg), s _ {a 2 s GG) 



2 6 7T 6 



da 



dp 



(1 - a - pf[{a + p)m 2 - 3aps}[(a + (3)m 2 c - aPsf 



d 4l m 



(l-a-(3) 2 m 2 c 



3q 



[2 (a + p)m 2 c - 3a/3s 



1-a-p 



[(a + P)m\ - 2a/3s][(a + (3)m 2 c - a/3s]} , 



P { i q) \s) 



(A4) 



f x ^(M 2 ) 



m 2 c (qg s a-Gq)(qq) f 1 da 2m 2 c - a(1 ^ )M s 



6tt 2 



a aM% 



n 



<GGG W , _ (9 3 s fGGG) [ x da f 



'(Mb) = - 



3 x 2 6 tt 6 J 



d0 



{ 



ln(a/3(l - a - /3)M| ) - 2 In 2 - -y E 
a/3 



[12(a/3A/|) 2 + 6a/3Mf (a + /3)m 2 c + {a + P) 2 m\ 



33(a(3M 2 ) 2 + 12a/3M|(a + /3)m 2 c + (a + P) 2 m\ (1 - a - Pfm\ 



np 



[2af3M 2 B + (a + /3) 



2a 3 



For 774: 



[3a(3M B + M 2 (a + P)m 2 c ]}e . 



P r\s) 



2 8 7T 6 



pf^) = 0. 



pf G \s) 



(9 2 S GG) 

2 8 7T 6 



da 



— r / -£(1 - a - /?) 2 [(a + /3)m^ - 3a/3s][(a + /3)m 2 - a/3* 



5/ 4tt 



(I - a- P) 2 m 2 c 



3q 



[2 (a + /^m 2 . - 30^5 



1-a-p 



[{a + f3)m 2 - 2a/3s][(a + P)m\ - a/3s]} , 



pT x {s) = 0, 



(A5) 



n 



mi x (qq)r M 2^ _ ^(gg^g ■ Gq){qq) f 1 da f m 



n 



(GGG) n , 2 



(M 



12tt 2 



(. 9s 3 /GGG) f 1 - f fi " 



l)e" 



3 x 2 9 tt 6 7 



da 



{ 



ln(a/3(l - a - /3)M|) - 2 In 2 - 7 _e 



[12(a/3M£) 2 + 6aPM 2 (a + p)m A c + {a + P) 2 m A} 



33(aPM 2 ) 2 + 12a/3Mf (a + P)m\ + {a + P) 2 m\ (1 - a - P) 2 m? c 



aP 



[2aPM% + {a + p) 



2a 3 



(o + /3)r 



[3a/3M| + M|(a + £)m 2 ]}e af)M B . 



2. The spectral densities of the interpolating currents with the quantum numbers J PC = + 
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For 775: 



P r\ S ) 

pf G \s) 



1 



2 7 7T 6 

0, 

(g 2 s GG) 

2 7 7T 6 



da 



" ^(1 -a- f3) 2 [(a + p)m 2 c ~ 3af3s][(a + [3)m 2 c - a/3s] 3 , 
f Pma * ,„ r {l-a- f3) 2 m 2 



70 m- 



3a 



-[2(a + /3)m 2 - 3a/3s 



l-a-/3 



[(a + (3)m 2 c ~ 2af3s][{a + f3)m 2 c - a/3s}} , 



m 2 (qq)' 



v/1 -4m 2 /s, 



(A6) 



n 



rnix(qq) . 



{M%) = 



n 



(GGG) 



(M 



m 2 (qg s a ■ Gq){qq) f 1 da , 2m 



12tt 2 



— (^T7T + 1 ) e 
a am ?> 



{ 



3 x 2 8 tt 6 J 
ln(a/3(l - a - /3)M|) - 2 In 2 - lE 
a/3 

33(a/3Af 2 ) 2 + 12a/3M 2 (a + (3 )m 2 c + (a + (3) 2 mj + (1 - a - /3) 2 m 2 c ^ + (a + ^ 



[12(a/3Af|) 2 + 6a/3Af|(a + /?)m 2 + (a + /3) 2 m^ 



a£ 



2a 3 



[3a/3A^ + Af^(a + /3)m 2 ]}e <^ M S . 



For r/ 6 : 

pf>( S ) 
pf G> ( S ) 



2 6 7T 6 

0, 

(g 2 GG) 

2 6 7T 6 



" ^(1 - a - /3) 2 [(a + /3)m 2 - 3a/3s][(a + /3)m 2 - a/3s] 3 , 



-^[2(a + /3)m 2 -3a/?4 



/Jmi„ 



(1 - a - P) 2 + 2a(3 
I/3 2 



[(a + (3)m 2 c - 2a(3s][{a + (3)m 2 c - a/3s] , 



2m 2 (qqY 



\J\- Am 2 c /i 



(A7) 



n™ x ^(M|) 



n 



(GGG), 



(M 2 B ) = 



m 2 (qg s a ■ Gq)(qq) f 1 da m 2 
Jo o?M 2 B 

3 x2V J J 



{ 



ln(a/3(l - a - fi)M%) - 2 In 2 - j E 
a/3 



[12(a/3A/ 2 f + 6a/3A/ 2 (a + l3)m 2 c + {a + /3) 2 m 4 c 



33(a/3A/ 2 ) 2 + 12aPM 2 {a + f3)m 2 c + (a + $) 2 m\ 3(1 - a - P) 2 rrt c 



a/3 



[2af3M 2 +{a + /3) 



P) -[3aj3M B + M 2 (a + /3)m 2 ] + s[3aj9Af| + M% (a + /3)m 2 ]}e 



2a 3 



1 - a - I3 l 
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For 777: 



P7 ert ( s ) = 1 f"^ f m "|(l- a ^) 2 [( a + ^- 3 a H[(a + /3K-^] 3 , 



2 5 7r 6 7 Q _ a 3 Jp m . n /? 3 
p<">(s) = 0, 

(GG), s (9 2 ,GG) f a — da r 0ma * o Al-a-P) 2 ml.. Q . 2 o fl . 

+ 15(1 ' Q 16 ^ 2 + 3 ° a/3 [(« + ^)m» - 2aM(a + " . 

pr x («) = 0, 

P r ] \ S ) = -^^-y^^Ts, (A8) 
mix{q q ) (M 2^ . 2mg(gff s g-Gg)(gg) da w 2 - a(1 j Mi 



JO 

ln(ag(l - a - flM£) - 21n2 - 7g [12(a/gM|)2 + 6 ^ m|(q + ^ + (q + ^ 
ap 

, 33(«/3AfI) 2 + UapM%{a + p)m\ + (a + pfm^ , (1 - a - /?) 2 m 2 + + 



a/3 

-[3a/3M^ + M|(a + /3)m 2 ]}e" 

For 77s : 



(1 ^L [^M% + M%{ 



Ps & = ^6 / I ^(1 - a - /3) 2 [(a + /3)m 2 - 3«M(a + - 



/3„ 



(*) = 0. 



^ ^ pja^ (!-a-^ 
3 x 2 b 7r b J Qm . n a^ Jp min a 

+ 3(1 ~" g/ f 2 )2 + 6a/3 [(« + /3)m 2 - 2a/3 S ][(a + p> 2 - a/3s}} , 

p™*(s) = 0, 

A* >"(-) = -"-^VT^US, (A9) 

Tf"te<to)(M*\ m 2 c (qg s a-Gq)(qq) f 1 da m 2 - Q(1 _"I% f . 

il§ (A7s) " 3^ Jo^K 



Mggfl - a - pOM|) - 21n2 ~ 7£ [12(a ^ )2 + 6 ^ m|(q + ^ + (q + ^ 
ap 

, 33(a/3M 2 ) 2 + 12a/3M| (a + /3)m 2 + (a + /3) 2 m^ (1 - a - /?) 2 mg ro 2 2 

1 -. |2apM B + (a + pjmj 



07 



(l-a-j8) 



[3a/3A/| + A/ 2 (a + /3)m 2 ]}e . 
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For r]g: 



pT{s) 



2 8 7T 6 



pp\s) = 0. 



pf G) (s) 



(g 2 s GG) 

2 8 7T 6 



J f &1 - a - /?) 2 [(a + ff)ml ~ 3a(3s][(a + f3)m 2 c - a(3s} 3 . 



(l - a- (3) 2 m 2 



3a 



■[2(a + 0)m 2 c -3a/3s 



l-a-/3 



[(a + /3)m 2 c - 2a(3s]{{a + (3)m 2 c - a(3s}} , 



P T x {s) = 0, 



pT> 00 



67T 2 



a/1 - 4m 2 /s , 



(A10) 



n 



mix(qq) ^yp. 



B) 



n 



(GGG) n , 2 



(M 



For ry 10 : 



m 2 c (qg s <7 -Gq)(qq) f l da m 2 c ~ a(1 ZW 



a aA/J, 



{ 



12tt 2 

3 x2V J 
ln(a/3(l - a - g)M|) - 2 In 2 - -yg 
a/3 



[12(a/3A/|) 2 + 6a/3Aff(a + /3)m 2 + (a + /3) 2 m^ 



33(a/3Af|) 2 + 12a/3M|(a + /3)m 2 + (a + P) 2 m 4 c (1 - a - /3) 2 ml 



a/3 



[2a0M% + {a + (3) 



2a 3 



[3a/3Af| + Af 2 (a + /3)m 2 ]}e "^"IT . 



2 5 7T 6 



da 
a 3 



(g 2 s GG) 

2 6 7T 6 



(1 - a - ^) 2 [(a + f3)m 2 c - 3a/3s][(a + /3)m 2 - a/3s] 3 , 



5/ 



2(1 -u~ P) 2 m 2 



■[2(a + (3)m 2 - 3a/3s] 



5(1 - a - /3) 2 + 10a/3 
I/3 2 

6(1- a -j3) 



[(a + /3)m 2 - 2a/3s][(a + /3)to 2 - a/3s] 



[(a + fi)m 2 c - 2a/3s][(a + /3)m 2 c - a/3s]} . 



P?o lx (s) = 0. 
P[l q) \s) : 



AmKqq} 2 



(All) 



n™* <9 " 9) (M 2 ) 



n 



(GGG) ( f./r1 



10 



2m 2 c (qg s a -Gq)(qq) f 1 da m 2 - n[1 _" c m 2 



o a aM 



(Mb) = - 



3 x 2 6 tt 6 



da 



d/3 



{ 



ln(a/3(l -a-^Ml) -21n2-7jj, 
a/3 



[12(a/3A/ 2 ) 2 + 6a/3A/ 2 (a + /3)m 2 + (a + ^) 2 m*] 



33(a/3A/ 2 ) 2 + 12a/?Mf (a + (3)m 2 + (a + fi) 2 m\ 3(1 - a - P) 2 m 2 c 



3(1 - a-/3) 2 
2^ 



a/9 

[3a/3M| + M|(a + /3)m 2 ] + 



[2af3M 2 +(a + /3) 



T~ 



1 - a - /3 



[3a/3Af^ + M| (a + /3)m 2 ] }e 



